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VO ■ Abstract 

O 

o, 

■rj" I We study electronic transport properties of disordered polymers in a quasi- 

- ^ ' one-dimensional model with fully three dimensional interaction potentials, 

"c^ ' We consider such quasi-one-dimensional lattices in the presence of both un- 

H . correlated and short-range correlated impurities. In our procedure, the actual 

— Jj . physical potential acting upon the electrons is replaced by a set of nonlocal 

Ph ', separable potentials, leading to a Schrodinger equation that is exactly solvable 

j^ I in the momentum representation. By choosing an appropriate potential with 

the same spectral structure as the physical one, we obtain a discrete set of al- 
gebraic equations that can be mapped onto a tight-binding-like equation. We 
then show that the reflection coefficient of a pair of impurities placed at neigh- 
ed I boring sites (dimer defect) vanishes for a particular resonant energy. When 
there is a finite number of such defects randomly distributed over the whole 
lattice, we find that the transmission coefficient is almost unity for states close 
to the resonant energy, and that those states present a very large localization 
length. Multifractal analysis techniques applied to very long systems demon- 
strate that these states are truly extended in the thermodynamic limit. These 
results are obtained with parameters taken from actual physical systems such 
as polyacetylene, and thus reinforce the possibility to verify experimentally 
theoretical predictions about absence of localization in quasi-one-dimensional 
disordered systems. 

PACS numbers: 73.20.Jc, 73.61.Ph, 73.20.Dx, 72.20.-i 
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I. INTRODUCTION 

Transport properties of disordered systems have become a fascinating research topic 
since the generahtv of locahzation phenomena in one dimension (ID) was first questioned 
a few years ago.Erta Opposite to the conventional view that in ID random systems almost 



all eigenstates are exponentially localized (see, e. g., Ref. |T9| and references therein), it is 
nowadays known that in disordered systems where disorder exhibits some kind of spatial 
correlation bands of extended states arise. Spatial correlation means that random variables 
are not independent within a given correlation length or, equivalently, that the noise is 
non-white. Furthermore, supression of localization by structural correlation has been found 
both in classical and quantum systems. In the quantum case, electronic transport has been 
of course the subject of most works. There exists at present much evidence that correlated 
disorder inhibits wave localization, and that bands of extended states appear in tight-binding 
Hamiltoniansai3'Ery as well as in more elaborated multiband systems like those described 
by Kronig-Penney models.ll3o Similarly, occurrence of super-diffusionia and reflectionless 
spin waves in Heisenberg chainal^ have been recently reported. In the classical case, random 
harmonic chains also present a band of shortwave delocalized vibrations whenever correlated 
disorder occurs,li3ll3 giving rise to a strong enhancement of the thermal conductivity of the 
lattice.li3 AH these theoretical analyses clearly demonstrate that transport properties in 
random systems where structural correlations are present are very different to what is found 
in ordinary random systems. It is also clear that supression of localization does not depend 
on the classical or quantum nature of the system and therefore structural correlations are 
to be regarded as the origin of this unexpected feature. 

In spite of the already available body of theoretical work, the physical relevance of these 
extended states is still unknown. To our knowledge, there is no experimental evidence what- 
soever of the existence of these states and their influence in measurable transport properties. 
We regard this as the key question to be posed about the theoretical results, above the more 
fundamental one on the nature of these states in infinite systems. Notwithstanding, we have 
addressed both issues along our research project on disordered systemsllHn!jOll3Ejc3 and 
in particular in the work we are reporting here. Regarding experimental demonstration of 
delocalization, we have recently shown how the bands of extended states must reveal them- 
selves through characteristic features in the dc conductance of disordered superlattices at 
finite temperature.tZl Moreover, we have also found that short-range correlated disorder has 
profound effects on coherentcd and incoherentc3 trapping, as well as on optical propertiesHll 
of excitons. In the same spirit, we have even proposed mechanical analogues where classical 
extended vibrations should be found.Ej We note that, aside from the basic research goal of 
finding out whether delocalization actually occurs in real physical systems, we also have in 
mind an applied aim, namely verifying whether correlated disorder gives rise to particular 
features that can be used for new devices or applications. 

In this work, we report further progress along the lines in the preceding paragraph. 
Searching for physically realizable systems where delocalization may play a crucial role, we 
turn ourselves to one of the pioneering works in the field, namely the work of Wu and Phillips 
on polyaniline (see Ref. | and references therein). These authors showed that polyaniline 
could be mapped onto a tight-binding, random dimer model that has a band of extended 
states, originated by a resonance of a single dimer defect. It is evident that on the one 



hand, similar mappings can be worked out for different polymers, and on the other hand, 
that delocalized bands have to affect their conductance properties. Indeed, Wu and Phillips 
argued that the fact that polyaniline was a conducting polymer was closely related to this 
unexpected delocalization phenomena.Q Their calculations were carried out in the framework 
of a purely ID tight-binding Hamiltonian. However, although it is quite reasonable to 
approximate the structure of a polymer by a line, it is also true that the physics involved 
is three-dimensional (3D) and that actually the linear structure of polymers is not straight 
but folds and wanders in 3D space. It is then natural to ask whether the above theoretical 
results will still hold when more realistic models including 3D effects are considered. The 
answer to this question is very important if Wu and Phillips's results are to be compared to 
measurements on real polymers: If the delocalized band is destroyed by 3D effects, then their 
theoretical results are merely academic, and worse, any possible technological application 
becomes very unlikely. On the other hand, this question is not without interest from the 
fundamental viewpoint. As we mentioned above, localization of almost all eigenstates by 
uncorrelated disorder is expected in ID random systems, but three-dimensional (3D) systems 
require a minimum amount of disorder to give rise to localization.llHl Then a question arises 
in a natural way within this context, namely the possible effects of correlated disorder on 
3D eigenstates. As far as we know, this problem has been already studied by Stephens and 
Skinner,B3 who found that tight-binding Hamiltonians with short-range correlated diagonal 
disorder in a cubic lattice present a localization threshold that is independent of the amount 
of correlation. This finding seems to indicate that the influence of structural correlations 
is relevant only to pure ID systems, which adds further interest to the elucidation of the 
applicability of Wu and Phillips's results. 

We address the above issues by introducing a new and completely general model to study 
electronic properties in random systems based on the so-called nonlocal (separable) potential 
(NLP) method, in which the actual potential at each site of an arbitrary lattice is replaced 
by a projective operator.c^Ha The treatment is fully 3D, although we restrict ourselves to 
a linear chain, so our model is not a ID model in the traditional sense as we have here an 
array of 3D potentials along a straight line. Moreover, the model can be straightforwardly 
extended to folded (i.e., non-straight) systems by appropriate choices of the parameters. 
As a major point, we will demonstrate the occurrence of a well-defined band of extended 
states in the electronic energy spectrum due to structural correlations, in spite of the 3D 
character of the equation of motion. In particular, we consider the case in which pairs of 
impurities (the so called dimer defects) are placed at random in an otherwise perfect lattice. 
The location of the band of the extended states in the electronic spectrum is determined 
from the condition of vanishing of the reflection coefficient from a single dimer defect. 

The paper is organized as follows. In Sec. II we present our model and summarize 
previous work of usB^ that is necessary for a better understanding of the present paper. 
After describing how nonlocal potentials can be used to model 3D systems, we consider the 
scattering from a single dimer defect exactly and find the resonance condition for perfect 
transparency. We close this section with a brief account of exact expressions to compute the 
physical magnitudes of interest in a lattice containing a certain number of dimer defects. 
Afterwards, in Sec. Ill we concern ourselves with our main topic, the random quasi-one- 
dimensional lattice with paired disorder. We present our numerical results demonstrating 
the existence of extended states via transmission and Lyapunov coefficients, which relate 



to physically relevant quantities such as localization length, as well as multifractal analysis, 
which points out the character of these states in the infinite size limit: Thus, we show how the 
bands of extended states reveal themselves through well defined peaks in the transmission 
coefficient versus energy plots, and how the truly extended character of those states is 
also demonstrated in their scaling properties. Section IV concludes the paper with a brief 
summary of the main results and some possible applications in different physical contexts. 

II. NLP APPROACH TO MULTICENTER INTERACTIONS 

The starting point for the NLP procedure is the Schrodinger equation for multicenter 
potentials, corresponding to the physical situation we want to gain insight into. The solution 
of this kind of problem is of widespread interest not only in condensed matter physics 
but also in atomic or molecular physics. As is well known, such solution is expected to 
involve enormous intricacies since in most cases prohibitively cumbersome calculations are 
needed. Several methods have been developed to study the motion of electrons in a given 
superposition of 3D potentials. Among them, the NLP approach is the natural generalization 
of the famous Kronig-Penney modelBZl to the 3D case. This method leads to an exactly 
solvable Schrodinger equation from which the electron energy can be obtained in a closed 
form. What is more important, it is always possible to find a NLP (or a sum of them) 
which reproduces any set of given electronic states,ca so there is no theoretical limit to the 
numerical accuracy with which physical results may be obtained. We first summarize the 
NLP formalism and then discuss its application to conducting polymers. 

A. Schrodinger equation for NLP 

We begin from the Schrodinger equation for NLP, which reads as foUowsEj (we take 
h = m = 1 hereafter) 

(p2 - 2E) ^(r) = Y: hV{\r - R,|) / rfVV(|r' - Rfc|)^(rO, (1) 

k ■' 

where R^ denotes the position of each lattice site and \k is the corresponding coupling 
constant. We will immediately see how Eq. (||) connects with the physical problem of 
interest, through suitable choices of the potential V . For simplicity we have assumed that 
the function V is spherically symmetric, although more complicated symmetries can be also 
easily handled. In Fourier space we have 

V'(P) = L2 ^ 2e ) ^ ^'^^^^ exp(-zp ■ R,)xfc, (2) 

where 

Xk = j d''pV\p)eM^V^k)Hv)- (3) 

Here '?/'(p) and V{p) denote the Fourier transforms of if) and V , respectively. The asterisk 
means complex conjugation; the Fourier transform of real and spherically symmetric func- 
tions is also real, but we retain it for if non-spherical functions are to be considered. The 



coefficients Xk are the quantities of interest, since we will show that they are related to the 
wave function in real space. We will be more specific about their meaning after we have 
specified the potential V{p) and computed the corresponding equations for Xk- Inserting 
Eq. (H) in Eq. (^ we obtain the following set of algebraic equations for the parameters Xk 



Xk = Yl \ j d'P 4^^ exp[^p • (Rfe - n,)]xr (4) 



Due to the spherical symmetry of the potential, the angular integration can be carried out 
in Eq. (^). In so doing, we finally obtain 

where Rkj = |Ra: — R-jl- It is understood that the factor {sin pRkj)/pRkj is replaced by 1 
when k = j, that is, by its limit as Rkj — > 0. 

B. Application to quasi-one-dimensional polymer models 

At this point we should stress that Eq. (^ is completely general, once the potential V{p) 
is specified. For a given ID, 2D or 3D lattice {Rfc, k = 1,2,... N}, N being the number 
of sites, the eigenenergies can be found by solving the secular equation arising from the 
N X N symmetric determinant associated to Eq. (^. In this fashion, we arrive at the key 
of the NLP procedure: The crucial question is to set up an appropriate potential V{p) 
that reproduces the observed energy values of the physical situation being considered. For 
instance, we have previously found that the Yamaguchi's NLPE3 is most appropriate to 
describe Coulomb bound states (see Ref. |26|), whereas surface ^-function potentials, that is, 
a force field vanishing everywhere except on a spherical shell of radius R, are very well suited 
to simulate electron potentials in long quasi-one- dimensional polymers, as polyacetylene or 
polyaniline.c3 This is the case we are interested in, and therefore we concentrate ourselves 
in this choice of potential from now on, i. e., we take 

V{r) = l,S(r-R), V(p) = ^'^. (6) 

Plugging this potential into Eq. (^ we get 

Xk = 8^ A{E)xk + E 8 4 B^mx,. (7) 

^ j^k ^ 

where for brevity we have defined 

sin^ pR n 



A{E)= dp^^^^J^ = ^(l-e-^^^) 



p2 - 2E Ak 

n (T?\ [^ A sin^ pR sinpRkj 7rexp{-KRkj) i^ d 1^ ^Q^ 



We are restricting ourselves to the case of interest, namely E < and then k = \J —2E is a 
real parameter. 

We note that interference effects due to the interaction of the electron with the lattice ap- 
pear in the coefficients Bkj{E): The larger the distance between site k and j, the smaller the 
corresponding coefficient. In other words, such coefficients are rapidly decreasing functions 
of Rkj whenever k is not very small (deep potentials). In our problem, this is a good approx- 
imation, and hence, to simplify numerical analysis, we assume that only nearest-neighbor 
interactions along the linear lattice are significant and write 

Xfc = 8 ^ A{E)xk + 8 ^ Bkk+i{E)xk+i + 8 ^ Bkk-i{E)xk-i. (9) 

To evaluate the coupling constant Afc in terms of experimentally measurable quantities, 
we consider an isolated potential centered at Rfc. This we accomplish by neglecting the 
interaction with other lattices sites so we take Bkk±i{E) ^ in Eq. (P), thus obtaining the 
condition determining the energy of bound states Ek, namely 8XkA{Ek) / R"^ = 1. This is 
a transcendental equation which can be easily solved numerically, but it leads to a simpler 
expression for small values of R, a limiting case we will consider later. For small radius R 
we can expand A{Ek) to obtain a relationship between the coupling constant A^ and the 
energy of the (single) bound state Ek = —i^l/2 

Afc~— (l + ztfei?). (10) 

Finally, inserting Eq. (^) in Eq. (|^) and taking the limit i? — ;> in such a way that Hk 
remains constant, we obtain the following tight-binding-like equation for the coefficients Xk 

exp{-KRkk+i) , exp(-Ki?fcfc-i) ^^^. 

[K-Kk)Xk = 5 Xfc+i H B ^fe-i- (11) 

Note that in this tight-binding set of equations the transfer integrals depend exponentially on 
the distance between nearest-neighbor distance, as it should be expected. This is consistent 
with our previous disregarding of longer range interactions. We can use these equations 
to describe the dynamics of electrons in the presence of diagonal as well as off-diagonal 
randomness. In the rest of the paper, and without loss of validity, we further assume 
diagonal disorder, which implies that Rkk±i = L, L being the lattice parameter. Defining 
p = kL and pk = HkL for the sake of brevity, we thus arrive at 

{p- pk)e''xk = Xkk+i + Xkk-i- (12) 

It is most important to stress that the number of free parameters appearing in these equations 
of motion has been kept to a minimum: We have only introduced the strength of the potential 
(which manifests itself in the value of the single bound state level, appearing in pk) and the 
lattice parameter L. 

Before we proceed to study the above equations of motion for our 3D model polymer, 
we now clarify the physical meaning of the coefficients Xk- From their definition in Eq. (^ 
and the Parseval identityl^ we have 



Xk 



J d\V{r)^{r + Kk), (13) 



with V{r) = {l/r'^)S{r - R). In the hmit i? ^ one gets V{r) -^ {l/r^)S{r) = S{r). 
Therefore, in that hmiting case 

Xk = ij{Tik). (14) 

We thus see that Xk is nothing but the value of the electron wave function at site R^, which is 
of course the quantity of interest. On the other hand, this is the reason of the denomination 
of Eq. ([l^ as tight-binding equations of motion. 

C. Perfect quasi-one-dimensional lattices 

Before considering random lattices, it is instructive to study the case of perfect lattices, 
that is, those lattices with Pk = Po- Since there exists translational symmetry, Bloch the- 
orem holds and we look for solutions of the form Xk = U exp{iQLk), Q being the crystal 
momentum and U a constant. Inserting this solution in the Eq. ([I^) we readily obtain the 
dispersion relation 



cosQL= —{p- Po). (15) 

Real values of Q, obtained by usual numerical methods, give us the electron energy E = 
— p^/2L^ as a function of Q and, consequently, the band structure of the lattice. To check 
the validity of the tight-binding approach we have assumed, it is necessary to compare this 
band structure with that obtained by including all non-nearest-neighbor interactions. We 
have already calculated it in Ref. ^, obtaining 

cosQ-^ = coshp — -e''°. (16) 

Assuming that the lattice parameter L is large and the potentials are deep (the basic as- 
sumptions in the tight-binding approach) it becomes clear that po and p are large but the 
difference p — po is small. With this assumptions it is a matter of simple algebra to demon- 



strate that Eq. (|T6D reduces to Eq. (^). This leads us to the conclusion that one can 
confidently use Eq. (|T5|) to describe the motion of tightly bound electrons in a lattice. 

From an experimental point of view, comparison with measurements of real polymers 
requires the evaluation of the two input parameters, namely L and po, from experimental 
data. The first one is usually known from X-ray data and, in principle, it is easy to obtain. 
The second one requires more information of the electronic band structure obtained, for 
instance, from spectroscopy measurements. Let us assume for the moment that the energies 
of the experimental band edges are known, and let Et and Eb be the energy of the top and 
the bottom of the band in the perfect lattice, respectively {QL = and QL = vr), and 



Pt = Ly/—2Et and pb = Ly/—2Eb. Using Eq. ([T5|) one has the following relationship 

Ptexp{pt)+pbexp{pb) , . 

Po = r^ r~^ (17) 

exp(pi) +exp{pb) 

Therefore, from the knowledge of the experimental band edges we can calculate the semi- 
empirical parameter po- We have thus shown how the model parameters can be found for 
comparison to the particular polymer one is interested in. We will make use of this result 
later. 



D. Scattering from a single dimer defect 

As mentioned in the Introduction, we are interested in the effects of structural correla- 
tions in the localization properties of quasi-one- dimensional polymer models. Following Wu 
and Phillips,u the simplest way to consider structural correlations is to introduce impurities 
at random but in pairs of sites. Physically this would correspond, for instance, to complexes 
of defects frequently encountered not only in polymers but also in molecular and solid state 
physics. In particular, a very clear description of dimer defects in polyaniline can be found 
in Ref. |^. In our model this means that p^ can take only two values, po and pg, with the 
additional constraint that pg appears only in pairs of neighboring sites, which we will refer 
to as dimer defects. 

Let us consider a single dimer defect placed at sites /c = and fc = 1 in an otherwise 
perfect lattice. To proceed, we have to take into account in the first place the condition for 
an electron to move in the perfect lattice which, recalling Eq. (|T3p, is given by 



■y (P - Po) 



< 1; (18) 



this constraint gives the allowed energy values once po is fixed. Now considering the equation 
of motion, Eq. (p!2D, at A; = —1, 0, 1 and eliminating xo and Xi o^i g^ts 

- X2 = (fi + fi' - fifi")X-i - (1 - fi")X-2, (19) 

where we have defined Vt = e''(p — po) and Vt' = e^{p — p^) for brevity. Besides a constant 
phase factor of vr, Eq. ([l^) reduces to the equation of motion in the perfect lattice whenever 
fi' = 0, in which sites k = and k = 1 have been eliminated. This means that the reflection 
coefficient at the single dimer vanishes, and consequently there exists a complete trans- 
parency. This occurs only for a particular energy of the incoming electron Er = — p^/2L^, 
given by the condition Q' = 0, i.e. pr = p'q- Hence this resonance effect occurs whenever the 
incoming electron matches the energy level of the (isolated) impurity, and this is possible 
only if Ipo — Po| < 2exp(— pg), as seen from Eq. (|^). This is to be compared with both the 



results of Ref. §, where a single resonant energy is found as well, and to those in Refs. ^ 



and IT3 for a continuum Kronig-Penney random dimer model, where an infinite number of 



resonances arise. We see that delocalization effects of structural correlation may be more or 
less dramatic depending on the physical situation studied. On the other hand, the important 
result is that the resonance of the simple ID tight-binding random dimer model is preserved 
in our 3D setup. 

E. Scattering from a lattice with random dimer defects 

We now proceed to the problem of a random lattice with a finite number of dimer defects. 
Of course, the above results do not imply anything about extended states in a lattice with 
a finite number of dimers defects, and it is necessary to study that problem separately. For 
definiteness, we introduce the concentration of defects c given by the ratio between the total 
number of impurities (twice the number of dimer defects) and the total number of sites 
N in the lattice. We introduce this definition to facilitate direct comparison with results 
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in ordinary random lattices with the same number of impurities and thus the same value 
of c, although in the latter case there are no constraints on the random location of the 
impurities. To study transmission properties of electrons through the random lattice, we 
place it between two semi-infinite perfect lattices. Therefore we introduce the reflection r 
and transmission t amplitudes through the relationships 

To determine both amplitudes we use the well-known transfer- matrix techniques (see, e. g., 
Ref. 0). Thus we cast Eq. ([I2|) into the matrix form 






(21) 
where a^ = (p — pk)e''. The transfer matrix of the whole lattice is then found as 

T{N) = n Pk, (22) 

k=N 

which relates the wave function at both edges of the lattice. Using the fact that det(T) = 1 
we finally arrive at the following expression for the transmission coefficient r = |tp 

, ^, 4sin^/9 , „. 

with 

D{E) = T^, + T^^ + Tl + r|2 + 2 (TnTi2 + T21T22 - T^T^i - T12T22) cosp 

- 2 (TnT22 + T12T21) cos^ p + 2 sin^ p, (24) 

where we have dropped the explicit dependence on N of the transfer-matrix elements. The 
transmission coefficient r can be recurssively computed from the matrix elements of T{N); 
taking into account the fact that T{N) = Pf^T{N — 1) and T(0) = Pq we find the following 
recurrence relations involving only real parameters 

Tn(iV) = aMT^,{N - 1) - Tn(iV - 2), 

Ti2(iV) = a^Ti2(iV - 1) - Ti2(iV - 2), 

T2i(iV)=Tn(iV-l), 

T22(iV) = ri2(iV-l), iV = 2,3--- (25) 

with the initial conditions Tjj(O) = 6ij, Tii(l) = ai, Ti2(l) = —1, 72i(l) = 1 and T22(l) = 0. 
Other physically relevant magnitudes can be readily obtained from the transfer-matrix 
T{N). In particular the Lyapunov coefficient, which represents the rate of the growth of the 
wave function, is nothing but the inverse of the localization length. It can be computed as 
(measured in units of L^^)E3 

nE) = (^) (Ti^i + T^2 + Tl + Ti,) . (26) 



Delocalization of electronic wave function is seen through the decrease of this parameter. 

The results we have obtained so far provide an exact, although nonclosed, analytical 
description of any random lattice with correlated as well as uncorrelated disorder. With 
them, we can compute the magnitudes we mentioned above. All expressions are very simple 
and suitable for an efficient numerical treatment for any specific cases. We will now evaluate 
them for several interesting cases to describe the relevant features of the model and the 
fingerprints of extended states. 

III. RESULTS AND DISCUSSIONS 

There are five parameters that can be varied in our model, namely the lattice constant 
L, the strengths of the scatterers po and p^, the total number of scatterers A^, and the defect 
concentration c. In order to find results as close as possible to actual systems, we consider 
a quasi-one- dimensional polymer, as is the case of polyacetilene (CH)^., which has been the 
focus of most of the experimental and theoretical works.E^ In the perfect lattice, taking an 
uniform carbon-carbon bond length of L = 1.39A,l23 we have previously estimatedc3 that 
Po = 1.466, corresponding to an energy level of the isolated potential of — 4.23eV. We can 
confidently take these values as correct since the predicted effective mass is found to be 
m* = 1.65, in excellent agreement with the experimental result m* = 1.7 ± 0.1. As an 
example, we will consider p^ = 1.550 implying an energy level of the isolated impurity of 
— 4.73eV. Note that the condition |po — Pol — 2exp(— Pq) holds, that is, the energy level of 
the isolated impurity hes in the band of the perfect lattice. Hence, according to our previous 
considerations, there exits complete transparency at an incoming energy of — 4.73eV if only 
a single dimer defect is placed in the lattice. Now we must elucidate what happens close 
to this resonant energy when several dimer defects are placed at random in the lattice, in 
comparison with lattices with the same number of unpaired defects. We used lattice sizes 
ranging from N = 2 000 to A^ = 500 000 sites. The largest of these systems are physically 
unrealizable, but it is important to study theoretically those systems to clearly elucidate the 
truly extended character of states close to the resonant energy, as we will demonstrate in the 
rest of the paper; the results for the smallest values are those directly related to experiments. 
Concerning the fraction of impurities, we only present here values corresponding to low 
defect concentration (c ranging from 0.1 to 0.3) because of their more physical relevance 
to actual systems, but we should stress that the main conclusion of the paper, namely the 
existence of truly extended states in quasi-one-dimensional lattices with correlated disorder, 
is independent of c. 

A. Transmission coefficient 

Since we are dealing with random lattices, we will need ensemble averages to compute 
the transmission coefficient. Some years ago, Sak and KramerE3 pointed out that only its 
logarithm obeys the central limit theorem, thus being the unique physically representative 
magnitude of the electron transmission, rather than the transmission coefficient itself or its 
inverse. Therefore we have actually computed exp (log r (£')), where (. . .) means ensemble 
average. Nevertheless, in what follows we refer this quantity simply as the transmission 
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coefficient and denote it by r, but it is understood that averages are carried out over the 
logarithm. 

An example of the behavior of r around the resonant energy is shown in Fig. |1], for both 
paired as well as unpaired lattices, with the same values of c = 0.1 and size N = 2000. 
A careful inspection of the figure clearly reveals that r is at least two orders of magnitude 
larger in paired lattices than in unpaired ones in the region of interest. In addition, and 
what is more apparent, r is close to unity around the resonant energy — 4.73eV, hence 
indicating that perfect transparency is preserved even when a finite number of dimer defects 
are placed at random in the lattice. This is a signature of the existence of a band of extended 
states close to that energy. We stress that, in spite of the fact that the plot corresponds to 
an average over 300 realizations, the transmission coefficient for typical realizations always 
behaves in the same manner, although plots are noisier. Thus, the only appreciable effect 
of averaging is to smooth out some particular very narrow, realization-dependent peaks, 
keeping the main common wide peak centered at the resonant energy. 

We want to highlight that the width of the transmission peak is always nonzero. Hence, 
close to the resonant energy, there is an interval of energies that also shows high trans- 
parency, similar to that of the resonant energy (note that there is a difference of about three 
orders of magnitude between the transmission coefficient of paired and unpaired lattices in 
that interval). The peak width depends on the concentration of dimers: The larger the con- 
centration, the narrower the peak, being always of finite width as already stated. Figure |^ 
shows results for two different values of c (0.1 and 0.3) and A^ = 2000. In addition, the 
width depends also on the system size: The larger the size, the narrower the peak, as shown 
in Fig. for two values of A^ (1000 and 2000) for c = 0.1. It is worth mentioning that r is 
always unity at the resonant energy, irrespective of the value of c or A^. 

B. Lyapunov coefficient 

The fact that the around the resonant energy r becomes close to unity suggests the 
possibility that the localization length of those states may be very large. This is, in fact, 
what is deduced from the analysis of the Lyapunov coefficient (recall that it is the inverse of 
the localization length). Results are plotted in Fig. ^ for paired as well as unpaired random 
lattices with A^ = 2000 and c = 0.1. The comparison between the results for the two kind 
of lattices is actually dramatic. First of all, we again observe that there exists a difference 
of several orders of magnitude between the localization length in both systems. In addition, 
paired lattices reflect the fact that a large number of states around the resonant energy 
present a very large localization length (which manifests itself in a deep minimum of F, with 
the Lyapunov coefficient taking values of the order of the inverse of the system size in a 
nonzero width region), whereas there is a monotonic dependence of the Lyapunov coefficient 
for unpaired disorder. 

C. Multifractal analysis 

From the study of the transmission coefficient and the Lyapunov coefficient we are led 
to the conclusion that there exists a large number of electronic states that remain unscat- 
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tered (or almost unscattered) by dimer defects. Such states are characterized by very large 
localization lengths (conversely, very small Lyapunov coefficients). Nevertheless, this result 
does not necessarily imply that those states are truly extended, namely states that cannot 
be normalized in the thermodynamic limit. Then we must search for a different approach 
in order to elucidate the localized or extended character of the eigenstates. The character- 
ization of the spatial extent of the wave function to all length scales may be accomplished 
by means of multifractal analysis.E^ Extended states span homogeneously the whole lattice 
whereas localized states remain confined in finite regions. The amplitude distribution of 
the electronic states can be characterized by the scaling with the system size of moments 
associated to the measure defined in the system by us (in our case the probability of finding 
the electron at a given point), 

yl^n=l \Xn\ J 

Notice that the second moment y^2(^) coincides with the inverse participation ratio (IPR) 
as introduced, for instance, in Ref. ^. The generalized dimensions Dg are determined via 
the scaling /ig(iV) ~ ]\[-('}-^)^i ^ for q ^ 1. For localized states Dg vanishes for all q whereas 
Dg equals unity for states spreading uniformly. In previous worksli3'll3 we have proven that 
multifractal analysis is a powerful tool to reveal the existence of truly extended states in 
ID random systems (phonons, electrons) with correlated disorder. Hence, we expect that 
similar characterization techniques also work well in quasi-one- dimensional systems, as in 
the present case. 

Let us start with the IPR. From its definition, it can be seen that delocalized states are 
expected to present small values of the IPR, of order of 1/A^, while localized states have 
much larger values. In the extreme case, when the electron is localized at a single site, 
Eq. ( p7D implies that fig{N) = 1. A typical situation is presented in Fig. ^ for the same 
system parameters as in Fig. |l|, using the initial conditions Xo = and xi = 1 to iterate 



the equation of motion given by Eq. (12) in order to find Xk- One can observe a wide, deep 
minimum of the IPR around the resonant energy for the paired disordered lattices, whereas 
this minimum is completely absent in the unpaired one. It is important to mention here 
that the value of the IPR at the minimum is independent of the defect concentration c, and 
depends only on the system size N. This suggest that the exact number of dimer defects is 
immaterial as far as the existence of such extended states is concerned. 

A complete multifractal analysis requires to study the scaling of all moments, defined 
by Eq. (p7|) with system size. We have consider such scaling for g = 2, 3 ... 6, and results 
are plotted in Fig. ^ for a concentration of c = 0.1 of dimer defects. We have observed that 
those moments scale very accurately as Hg{N) ~ A^~('?~i) for energies close to the resonant 
one, as illustrated in Fig. ^ for — 4.70eV (close but not exact the resonant energy). On the 
contrary, for more distant energies fJ^giN) follows a power law for small systems but tends 
to a constant value for larger ones, as seen in Fig. § for — 4.50eV. Therefore, according 
to the above discussion, the generalized dimensions Dg are exactly unity, within numerical 
accuracy, for states close to the resonant energy, thus indicating the truly extended character 
of such states, in agreement with results obtained from the analysis of the IPR above. 
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IV. CONCLUSIONS 

In this paper we have considered electron dynamics in quasi-one-dimensional models of 
polymers with correlated disorder, and we have compared our results to those obtained in 
systems with ordinary (uncorrelated) disorder. Our procedure based on NLP allows us to 
carry out a fully three dimensional analysis of the model, with the scatterers placed along 
a straight line. It is important to realize that this technique can be made exact so there 
are no theoretical limitations in this approach. In addition, the exact solution can be found 
for an arbitrary NLP, as we actually demonstrated [see Eq. (^)] by means of the Fourier 
transform (which, in fact, is completely equivalent to use a Green's function formalism). As 
our selection for a suitable NLP that can reproduce experimental data for polymers, we have 
used surface (5-function interactions with vanishing radius since, as we previously shown,ca 
this potential gives very accurate results in the context of polymers. Using this model, we 
have found that there exists a resonant energy for which the reflection coefficient of a single 
dimer defect vanishes, that is, there is perfect transparency. Afterwards, we turned to the 
problem of electron scattering when several of such defects are located at random in the 
lattice. Results from the evaluation of the transmission coefficient and Lyapunov coefficient 
(the inverse of the localization length) strongly suggest that there exist many states close to 
the resonant energy that remain unscattered, where this is not the case when the constraint 
of pairing is relaxed. To demonstrate that such states are actually extended in nature, we 
have used multifractal analysis, which confirms our claim. 

We now stress the physical relevance of our results. A key observation is that the resonant 
energy value does no depend on the defect concentration c. Therefore, by modifying this 
concentration, we could shift the Fermi level of the quasi-one-dimensional lattice to match 
the resonance. In this case, when the Fermi level reaches the resonant energy, a large 
increase should be observed in the electrical conductance peak. In fact, we have recently 
demonstrated in ID Kronig-Penney models with correlated disorder that very noticeable 
peaks in the finite temperature dc conductance appear as the Fermi level is moved through 
the band of extended states.O In a similar way, one could expect such dramatic increase 
in more elaborated models, as it is the case with the one we present here. On the other 
hand, our model supports the results previously found for simpler, pure ID ones.l3 It seems 
to us that this agreement makes very appealing the idea of use polymeric systems to confirm 
experimentally the existence of delocalized states in ID model, for what we have shown is 
that 3D effects do not destroy the coherence required for those states to appear. In fact, it 
is very tempting to relate all this to the known fact that polyaniline shows a metal-insulator 
transition, with the concentration of dopant acting as a control parameter (see Ref. ^ and 
references therein). Another important consequence of this work is that other experimental 
procedures we have proposed to find out whether delocalization can be measured or not, 
such as disordered superlattices,t3 are likely to be correct even if 3D effects have not been 
taken into account. We hope that this result encourages experimental work in the field of 
quasi-one-dimensional disordered systems which could give status of physically relevant to 
the theoretically predicted bands of delocalized states. 
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FIGURES 

FIG. 1. Transmission coefficient as a function of the energy around the resonant energy 
— 4.73eV, for paired (upper curve) and unpaired (lower curve) random lattice. Shown are av- 
erages over 300 realizations. Every realization consists of A^ = 2000 scatterers and a fraction of 
defect c = 0.1. 

FIG. 2. Transmission coefficient as a function of the energy around the resonant energy 
-4.73 eV for paired lattices with (a) N = 1000 and c = 0.1, (b) N = 2000 and c = 0.1, and 
(c) N = 2000 and c = 0.3. Shown are averages over 300 realizations. 

FIG. 3. Lyapunov coefficient as a function of the energy around the resonant energy — 4.73eV, 
for paired (lower curve) and unpaired (upper curve) random lattice. Shown are averages over 300 
realizations. Every realization consists of iV = 2000 scatterers and a fraction of defects c = 0.1. 

FIG. 4. IPR (^2) as a function of the energy around the resonant energy — 4.73eV, for paired 
(lower curve) and unpaired (upper curve) random lattice, with the same parameters as in Fig. 1. 

FIG. 5. Scaling of moments ^2 to ^e with the system size in paired random lattices for an energy 
— 4.70eV (solid lines), i. e., close to the resonant energy, and for an energy — 4.50eV (dashed lines), 
i. e., far from the resonant energy. Defect concentration is c = 0.1. 
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